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Abstract 

we explore the consequences of assuming a simple 3-parameter form, first 
without T-violation, for the neutrino mass matrix M in the basis z/g, u^, i/^ with 
a new symmetry. This matrix determines the three neutrino masses mi ,m2 ,m^, 
as well as the mapping matrix U that diagonalizes M. Since U, without T- 
violation, yields three measurable parameters Si2, S23, S13, our form expresses six 
measurable quantities in terms of three parameters, with results in agreement with 
the experimental data. More precise measurements can give stringent tests of the 
model as well as determining the values of its three parameters. An extension 
incorporating T-violation is also discussed. 
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1. Neutrino Mapping Matrix without T-Violation 

In this paper we wish to explore further the connection between the neutrino 
mass operator Ai which contains three neutrino nnasses mi, 7712, and the 
neutrino mapping matrix U, characterized by the standard four parameters 
^12, ^23, 013 and e^^. For clarity, we first examine the special case that the 
T-violating phase parameter 6 = 0. In terms of the mass eigenstates z^i, 1^2 
and 1/3 the neutrino mass operator is 

M = miViUi + 7712^2/^2 + rnsVsi^s- (1-1) 

Our assumption is that the same Ai, when expressed in terms of z/g, i^^ and 
Vr, has a simple form with a new symmetry property: 

also with three real parameters a, (3 and itiq. These three new parameters are 
to be determined by the mass eigenvalues mi, m2 and m,3. The transformation 
matrix U that brings A4 from (1.2) to (1.1) is the neutrino mapping matrix 
for 5 = 0. (The general case when 6^0 will be discussed in the next section.) 
Throughout the paper, we denote 

Ui = ij{iyi) and 77^ = ^t(^^.)^^ (^1_3) 

with ip{iyi) a 4-component Dirac field operator, ^ denoting the hermitian con- 
jugation and the index i = 1, 2, 3 or e, /i, r. 

Since the neutrino mapping matrix U is independent of the overall mass- 
shift term mo, in order for our hypothesis to be successful, there must be some 
special features about the first two terms in (1.2): 

(^{T^t - T^,,){yT - ^ii) + (3{y^,- Ve){p^ - l^e)- (1-4) 

We note that (1.4) is invariant under the transformation 

i/e^ i^e + z, i^f^^ i^f^ + z and Ur ^ I^t + z (1.5) 

with z a space-time independent constant element of the Grassmann algebra, 
anticommuting with the neutrino field operators z/j. Thus, the usual equal-time 
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anticommutation relations between the neutrino fields Vi and their zero-nnass 
free particle action-integral are invariant under (1.5). This synnnnetry is violated 
by the last mo-dependent ternn in (1.2), as well as by T-violation, as we shall 
discuss later. The interesting case that z might be space-time dependent will 
not be discussed in this paper. 

Expression (1.4) can be generalized to an equivalent form with three real 
parameters a, b and c: 

aiVr - V^){Ur - V^) + h[Vf, - Ve)iUi_i - Ve) + c[Ve - Vr)[Ve " ^^r)- (1-6) 

The corresponding neutrino mass operator is 

i 

It is clear that (1.6) is also invariant under the transformation (1.5). The same 
invariance can also be expressed in terms of the transformation between the 
constants a, 6 and c, with 

a— ^a + A, 6^6 + A, and c^c + A. (1.8) 

As we shall prove, the form of the neutrino mapping matrix U remains un- 
changed under the transformation (1.8). 

Since the relative phases between z/e, and Ut are unphysical, we may 
transform 

Ue -Ve, i^fi -1^1^ and z/^ i/r, (1.9) 
so that (1.7) is written in a less symmetric form, with 

(1.10) 

The sole purpose of using this less symmetric expression of M. is to have the 
resulting neutrino mapping matrix U in the standard form given by the particle 
data group[l]. We write (1.10) as 

_ ( uA 

M = (z7g 17^ 17^) (mo + M) 



V J 



(1.11) 



where 



M = 



I b + c -b c 
—b a + b a 
\ c a c + a 



(1.12) 



The neutrino mapping nnatrix U is defined by 

/ mi \ 



U\mQ + M)U = 



Introduce a 3 x 1 column matrix 



4>2 — 

One can readily verify that 



1 
^3 



m2 

V 7723 y 



/ 1\ 

1 

-1 



(1.13) 



M02 = 0; 

i.e., 02 is an eigenvector of M with eigenvalue 0. Let 
two real normalized eigenvectors of M. Since 

with ~ denoting the transpose, the neutrino mapping matrix U is 

U = {(1)1 (/»2 h), 

which, on account of (1.14) and (1.16), is given by 



(1.14) 
(1.15) 

)i and 03 be the other 



(1.16) 



(1.17) 



U = 



( 



\ 



cos I 



sin2 



gSin| + ^2C0s| 



yjcosf + yfsinf -y^ -y^sin f + cos f 



(1.18) 



in the approximation that the T-violating parameter ^ = 0, with the angle 
denoting the azimuthal orientation of 0i, 03 around the fixed eigenvector 02- 
Except for minor notational differences, the above U is the same expression 
first obtained by Harrison and Scott[2]. 



Next we return to the transformation (1.8), under which M of (1.12) trans- 
fornns as 



f 2-1 1 \ 
-1 2 1 
1 1 2 



Since 



/ 2-1 l\ 
1 2 1 
1 1 2 



V 



= 0, 



the neutrino mapping matrix U remains given by (1.18). Setting 

A = — c, 



(1.19) 



we have 



a 



a 



a 



c, 



I3 = b-c, 
0. 



(1.20) 



The corresponding neutrino mass operator A4 of (1.7) becomes (1.2). With 
the additional phase convention (1.9), M of (1.10) reduces to 



1.21) 



which has only three parameters a, (5 and mo. Of course, the mass operator 
(1.21) is a special case of the mass operator (1.10), which has 4 parameters 
a, 6, c and mo. It is of interest that they shares the same neutrino mapping 
matrix U given by (1.18), provided that a — c = a and b — c = (3. Yet, the 
neutrino masses mi, m2 and ms in the two cases can be different, as can be 
readily seen by examining the trace of M given by (112). Therefore, the full 
physical contents of (1.21) and (1.10) are not the same. This is especially 
important when we generalize the model to include T-violation in the next 
section. 
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For the remaining part of this section, we shall explore further the physical 
consequences of our nnodel, using only the more restrictive form (1.21) with 
three real parameters a, (3 and ttiq. 

It is instructive to re-derive (1.18) in a more elementary way. Write (1.21) 

as 

( vA 

M = (z7e Z7^ Vr) {aM^ + (5M^ + mo) 



(1.22) 



with 



/ \ 
1 1 
1 1 



(1.23) 



and 



(I -1 0\ 
-1 1 

^ ^ 



(1.24) 



(1.25) 



The matrix aMo: + jSMp in (1.22) will be diagonalized in two steps. Introduce 
first a real orthogonal matrix [3,4] Uq by setting 6* = in (1.18); i.e.. 

The matrix Uq diagonalizes Mq,, with 

/ \ 

= uImJJq = 2 



and transforms M/j to 




\0 1) 



(30 


-Vs 1 



(1.26) 



(1.27) 



Their sum aM'^ + /3Mp can then be readily diagonalized with another real 
orthogonal transformation matrix 



/ cos f - sin f \ 





sin I cos 1 



1 


2 



2 J 



with 



sin^ 
cos = 



and therefore 



1 — i 

tan^ = 



2a- p 

The resulting transformation matrix U = UqUi satisfies 



(1.28) 

(1.29) 
(1.30) 

(1.31) 



U 



(1.32) 



and is given by (1.18). The corresponding masses mi, and 777,3 ai'e related 
to a, (5 and 7770 by 

3/?2 



and 



?77i = a + - (a - -) 



7772 = 7770 



7773 = a + ;5 + (a - ^) 



1 + 



1 + 



(2a - ;5)2 



3/?2 



+ "7o, 



+ 7770. 



(1.33) 

(1.34) 
(1.35) 



(2a-/?)2- 

The matrix t/ depends only on one parameter Q, which in turn is determined 
by the ratio 

In the standard parametric representation, the matrix element [/13 is S13 = 
sin ^13 when e^^ = 1, with the experimental bound[l] 



4 = 0.9 g X lO-l 



(1.36) 
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From (1.18), Ui^ is — y|sin|. It follows then 



sm2^ = ^4«l- (1-37) 



Thus, by using (1.29)-(1.31) we see that 



q «1, (1.38) 



a 

which together with (1.33)-(1.35) yield the conclusion that mi and m2 are 
very close, forming a doublet, and ms is the singlet. Their mass differences 
are given by approximate expressions: 

m2-mi = -ll3 + 0(—) (1.39) 

ms-m2 = 2a + l /3 + o(—) (1.40) 

and ^ 

ms - hmi + m2) = 2a-] f3 + 0(—). (1.41) 

From mi < 1712, we conclude 

P <0. (1.42) 

Furthermore, z/3 is heavier or lighter than the doublet ui and 1^2 depending on 
the sign of a, with 

a > for 7713 > "T-i 01' 
and (1.43) 

a < for 7713 < "^1 01' 

Neglecting 0{f3/a) corrections, we have from (1.34), (1.39) and mi positive, 

mo>l m (1.44) 

and 

Sm'^ = ml-ml = (mo-- \/3\^3 \I3\. (1.45) 



Thus 

5m^ > ^ 0". (1.46) 



For 



Am^ = - i + ml), (1.47) 



we find, neglecting 0(/?^), 

Am^ = 4a (a + mo) + mo - 2a) \(3\. (1.48) 

The experimental values for dm?' and Am^ are given by[l] 

(^m^ = 7.92 (1 ± 0.09) x 10"W (1.49) 

and 

lAm^l = 2.4 (l ) X 10-=^ evl (1.50) 



Their ratio is 



= 3-3 1 X 10 ^ (1.51) 



|Am2| ■ V -0.28 

Next, we analyze first the case that the singlet 1/3 is of a lower mass than 
the doublet masses; i.e., a < 0. In that case, since ma > 0, (1-26) yields 

ms = mo - 2 |a| - i \/3\ + o(^) > 0; 

therefore 

mo > 2 \a\. (1.52) 
Neglecting 0{(5/a) corrections in (1.45) and (1.48), we have 

6m? 3 mo 



Am^ 4 a mo — la I 



(1.53) 



which gives 

3.(3. .5m\ 3. p. ,^ 
2 'a' 'Am2' 4 'a' ^ ^ 
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Combining this expression with (1.51), we find 

4.4 X 10~2 > 1-1 > 2.2 X 10~^ 
a 

On the other hand, fronn (1.29) and to the same accuracy, we have 

3/?2 



(1.55) 



sin^^ = 



4^2 ' 



(1.56) 



which on account of (1.36) gives 



a' 



^ = (0.72 



+1.84 
-0.72 



) X 10"^ 



(1.57) 



While (1.55) is barely consistent with (1.57), the compatibility depends on that, 
within one standard of deviation, (1.57) is also consistent with = (i.e., 

= 0). Thus, this "compatibility" between (151) and (1.57) is definitely 
not a comfortable one. A more accurate determination of Ui^ may well rule 
out the case that us can be lighter than the doublet z/i, z/2. Within our model, 
we also made a similar analysis for the case that the singlet is heavier than 
the doublet vi, v^. in that case, a > and the situation is quite different; 
there is no incompatibility between (1.51) and (1.57). 
Remark . We note that if /3 = in (1.21) then there is only one term 



(1.58) 



that is relevant for the determination of the mapping matrix; correspondingly, 
in the mass operator (1.22) we need only to consider aM^, with given by 
(1.23). Introducing a 45^ rotation matrix 



we have 



(1.59) 





( 1 


\ 


Ri = 


^ 






lo - 








^ \ 


RlMaRl = 









^ 2 ^ 



(1.60) 
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Because of the degeneracy in its first two eigenvalues, RiMaRi comnnutes 
with any unitary matrix of the form 



V 



u 



\ 


1 y 



(1.61) 



where w is a 2 x 2 unitary matrix. Thus there is a one-parameter family of 
solutions for the neutrino mass eigenstates. 
The situation is quite different when 



a 



= + 



(1.62) 



As mentioned before, because of the invariance (1.5) and the phase convention 
(1-9). 

^2= \^ -{l^e + l^r) (1-63) 

is a mass eigenstate. Furthermore, the transformation matrix 

Uo = R1R2 (1.64) 

is completely determined, with 



R2 = 



I 



\ 



1 



(1.65) 



which is a rotation of angle= sin ^ y^. For fija small but nonzero, the 
mapping matrix JJ deviates from ?7o through the small parameter as given 
by (1.18). 
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2. Neutrino Mapping Matrix with T-Violation 

We generalize the neutrino nnass operator Ai by inserting phase factors e^*^ 
into (1.6), replacing it by 

a(T7r - V^) {yr - z/^) + hiy^ - Ve) (i^/i - Ue) + c{e'^'^Ve - Vr) (e'^'z/g - ^t) (2.1) 

where a, 6, c and r] are all real. When = 0, (2.1) becomes (1.6), and is 
invariant under the symmetry (1.5). Furthermore, if e*^ ^ ±1, T-invariance 
is also violated. As in (1.6), in order to conform to the standard form of the 
neutrino mapping matrix U given by the particle data group[l], we make the 
phase transformation —Ue, —v^ and Vr ^t, the mass operator 

(1.10) is then replaced by 

c(e-*^Fe + F^)(e*Ve + l^r) + E^i^i, (2-2) 



which can be written as 




(2.3) 



where 



M = aMa + bMb + cMc + tuq 



(2.4) 



with 



a 



/ ^ 
1 1 



(2.5) 



\0 1 1) 



( 1 -1 o\ 
Mb= -1 1 
^ ^ 



(2.6) 



identical to Ma and Mp given by (1.23) and (1.24), and 

/ 1 e-"' \ 
Mc= 

. e^^ 1 , 



(2.7) 
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As in (1.25)-(1.27), we first perform the Uq transfornnation. Let 



1^,(1 + e-'^ - 2e'^) |(l-cosr7) ^(-1 + 



Next, we apply the Ui transfornnation given by (1.28), and write 

UiUoMUoUi = Ho + ch 
where Hq is diagonal, given by 

/ //I \ 



(2.8) 



(2.9) 



fl2 

\ fis ) 



(2.10) 



with /ii, /i2, the sanne ones in (1.33)-(1.35), except for the replacement of 
a, (3 by a, 6; i.e.. 



and 



IJLi = a + h - {a - 1 + 



112 = mo 



lis = a + b + {a - -) 1 + 



3b' 



(2a - 5)2 



36^ 
(2a - 6)2 



(2.11) 



+ mo. 



In (2.9) 

h = UiM'jJi. 

Since Uo and [/i are real and symmetric, /i is a hermitian. 
It is useful to decompose h into real and imaginary parts: 

h = h^ + ih^ 



(2.12) 



(2.13) 



13 



where 



= sin 77 



^cosf + yi 

-^cosf -V^sinf 



and the matrix elements of are given by 



^11-3 







(2.14) 



2 + - cos ^ + (1 + cos 0) cos 77 + \ Q ( o + 77 ) sin ^ 



^22 = ^(1 - cos 77), 



^33 = "^(2 + COS 77) — ^(1 + 2 COS 77) cos^ — ^ 



1 



6 



2^3 



-(1 + 2 cos 77) sin^, 



1 



^12 - ^21 - 3^ 2 



-(cos ^ — v^sin ^)(1 — cos 77), 



(2.15) 



/ifg = /i£ = i (^/3 COS ^ - sin ^) (1 + 2 cos 77) 
6 



and 



uR — un _ 
'hs — 'h2 — 



\ 7:(cos-H ^ sm -)( 1 — COS77). 

^6^ 2 V3 ^ 



The presence of violates T-invariance. 
We note from (2.14) that the element 



1 



\3 



sin77 



(2.16) 



is of particular importance for testing T-invariance. Furthermore, there are at 
least three cases to be considered: 

i) |c| << then T-violation is much smaller than the present upper limit, 
regardless of 77. 

ii) |c| ~ 0[|6|] but |sinr/| << 1, then T-violation is again very small on 
account of the prefactor sin 77 in (2.14). 
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iii) |c| ~ 0[\bW and | sin?7| ~ 0[1\; then T-violation can be close to the present 
upper linnit. 

The diagonalization of the 3x3 matrix (2.9) is simplified in case i). In that 
case, |c| is much less than |5| and \a\. The mass eigenstates and the correction 
to the neutrino mapping matrix can be readily obtained by using the standard 
first order perturbation formula. 

Another simple case is I77I << 1, which includes the above case ii). De- 
compose (2.7) into a sum 



Me = (Me)o + A 



(2.17) 



with 



(Me)o = 



(1 1\ 


1,1 ij 



(2.18) 



and 



A = 



e-^^ - 1 \ 

\e'i-l 



(2.19) 



Correspondingly, (2.4) can be written as 

M = Mo + cA 

with 



(2.20) 



Mo = aMa + bMb + c(Me)o + mo. (2.21) 

Mo can be diagonalized by the same unitary matrix (1.18), with the angle 6 
given by (1.29)-(1.31), in which a and /3 are given by (1.20). For \r]\ « 1, A 
is small; the neutrino mapping matrix U can then be derived by using (2.20) 
and treating cA as a small perturbation. 

We wish to thank W. Q. Zhao for helpful assistance and for informing us 
of the pioneering papers of Refs.[2,3], after our completion of this manuscript 
except for its references.. 



15 



References 

[1] S. Eidelman et al. (Particle Data Group), Phys. Lett. B592, 1(2004). 

[2] P. F. Harrison and W. G. Scott, Phys. Lett. B535, 163(2002). 

[3] L. Wolfenstein, Phys. Rev. D18, 958(1978); 

P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B530, 167(2002); 

Z. Z. Xing, Phys. Lett. B533, 85(2002); 

X. G. He and A. Zee, Phys. Lett. B560, 87(2003). 

[4] T. D. Lee, American Physical Society Meeting, First Session on 50 Years Since the 
Discovery of Parity Nonconservation in the Weak Interaction I, April 22, 2006, 
Chinese Physics 15, 1125(2006). 



16 



